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The charged particle environment in polar orbit can be of sufficient in- 
tensity to cause spacecraft charging. In order to gain a quantatative under- 
standing of Such effects, the Air Force is developing POLAR, a computer code 
which simulates in three dimensions the electrical interaction of large space 
vehicles with the polar ionospheric plasma.. It models the physical processes 
of wake generation, ambient ion collection, precipitating auroral electron 
fluxes, and surface interactions, including secondary electron generation and 
backscattering, which lead to vehicle charging. These processes may be fol- 
lowed dynamically on a subsecond timescale so that the rapid passage through 
intense auroral arcs can be simulated. POLAR models the ambient plasma as iso- 
tropic Maxwellian electrons and ions (0 + , H + ), and allows for simultaneous 
precipitation of power-law, energetic Maxwellian, and accelerated Gaussian dis- 
tributions of electrons. Magnetic field effects will be-modeled in POLAR but 
are currently ignored. 

The theoretical models and approximations employed in POLAR are dis- 
cussed, including an effective process for stabilizing the Poisson-Vlasov it- 
eration process in the short Debye length extreme. A preliminary POLAR cal- 
culation is presented which predicts the effects of measured auroral fluxes on 
the shuttle orbiter, and demonstrates the combination of conditions required 
for substantial differential charging of the orbiter. 

INTRODUCTION 

The charging of the space shuttle orbiter in the polar (auroral) iono- 
sphere has been investigated by Katz and Parks (ref. 1). In that paper, the 
authors argue that observed precipitating electron fluxes can exceed the ram 
ion flux (~1?0 pA/nr for n = 10 s cm, V = 8 km/ sec). This implies a 

possible overall current balance at a negative potential to enhance ion col- 

lection and retard electron fluxes. Both theoretical and experimental (refs. 

2, 3, 4) studies have shown that as object dimensions become large with re- 
spect to the local Debye length, space charge effects will severely reduce the 

ion current collection compared to orbit limited theory. This leads to a 
"size-effect" for the potential buildup on objects subjected to auroral elec- 
tron fluxes. Katz and Parks (ref. 1) have calculated, for a conducting 
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it necessary to include wave effects and turbulent phenomena by the use of 
time averaged models, diffusion rates, heating, etc. 


the lower auroral ionospheric environment is summarized in Table 1. The 
values given in the table are considered, typical, but none are fixed in POLAR 
where these and other parameters (such as the oxygen to hydrogen ion ratio) 
are variable over wide ranges. The energetic electron spectrum is currently 
modeled as a sum of power law, hot Maxwellian,, and Gaussian distributions.- 


POLAR presently ignores the magnetic field; however, its inclusion is a 
current effort. Some of the methods under study for including magnetic ef- 

TGCtS 


Place vxB potential gradient on the vehicle. 

Include B in ion trajectories interior to sheaths. 

Pitch angle conics for the energetic electrons (presently, isotropy is 
assumed). 

Modification of secondary and photoelectron emission where B parallels 
surface. H 

Modification of secondary an^. photoelectron surface conductivities in 
directions perpendicular to B. 


a 


We now outline the methods used in POLAR to perform charging calcula- 
tions. A calculation is broken in to the major steps listed below: 

Vehicle definition. 

Environment specification and computational grid construction. 
Presheath and wake ion density calculation. 

Initial surface charging using flux estimates. 

* Poisson and electron charge density calculation. 

* Sheath determination and particle (ions) tracking to determine: 

sheath ion densities 
ion surface currents 

* Surface charging. 

After all of the above modules have been executed once, the * items are 
iterated upon to produce a final solution. At this Stage in the development 
of POLAR, only negative surface potentials are allowed. Although slightly 
positive potentials could occur under natural charging conditions due to sec- 
ondary or photoelectron emission, the ambient electrons should limit positive 
potentials to a few kT at most. With this constraint, electrons are consid- 
ered to be repelled with densities given by the Boltzmann expression, n e = 
n o exp(-qV/kT). The methods used to calculate the attracted ion densities 
are described later in this paper. 


Positive potentials may also be achieved by the emission of electrons 
(ref. 6). The physics of electron collection in the lower ionosphere appears 
more complex than the collection of ions due to turbulent processes in both 
the emitted beams and in the ambient plasma (refs. 7, 8). Thus, spacecraft 
generated high positive potentials are not now considered in POLAR but will be 
addressed in future work. 
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This "neutral ion" approximation to the density is particularly^sirople 
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Pitaevskii (ref. 10) of the flow of a hypersonic plasma over a semi-infinite 

neutral ion approximation densities calculateo by POLAR- 

In fiaure 6 we present a comparison of the corrected neutral ion den- 

sity, calculated by POLAR, with electron density me WpTwgmuir 
et al (ref. 11). These measurements were made on STS-3 by the PDR cang 
Drobe* For STS-3, the estimated range of Mach numbers was 5-8, and the est 
mated* uncertainty in the absolute scale of electron densities was 2-5 mM 
densities were calculated at Mack numbers of 6 and 7 using the shuttle moae 
Dresented in this paper.- A comparison between electron and ion densities is 
ctrirtiv valid onlv when the plasma is quasineutral and the actual degree to 
whkh the wake plasma is quasineutral cannot be determined without a complete 
analysis* however the deviations from quasineutrality are probably within the 
experimental error and numerical uncertainties. The emph as is °f th^s eompar - 
son is placed upon the horizontal agreement which demonstrates POLAR s ability 
to model the wake edge using actual orientation data. 

As a calculation proceeds, POLAR locates (by inspection of potentials) 
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placement process is currently being implemented in POLAR, the calculations 
presented here used only the neutral ion densities. 


SURFACE CHARGING 
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no . n models the accumulation and transport, of charge on the vehicle 
g umped element circuit analogy. In this equivalent circuit each sur 
face cell ana conductor represents a node of the circuit. Surface voltaoe« Ur * 

eqSation t6d ^ imp11citly tim estepping a differenced approximation to the 
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f? d ~ are current and voltage vectors with each surface and conductor 
rK^r,ecto7rlp§s5 2* Capacita "“ 2 »• conductance. 
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where we have the current due to ions, ion impact secondary electrons auroral 

electrons* ^hTauror^T ? econdary and backscatter electrons* and photo- 
etectronsi The auroral electron currents and their secondary and backscattpr 
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near’aero until U d aia<S i^! U n e °f. the, r o”ission and the surface is held 
zero until it again shows negative charging behavior. This techniuue 

prevents the oscillations that could occur with a combination of larqe time 
step and an extremely voltage sensitive current source. 9 


As mentioned previously, POLAR models the flux,* of the eneroPtir 

d^s tr i but i ons^ rpf 3S 1 ? \ C0 ^ * na ?; ’ ’’I!/ p ? wer law » hot Maxwellian and Gaussian 
distributions (ref. 12) given by the following expression: 
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where K is kinetic energy, V is the surface potential, F « n •/kT/2itm, and A 

hasten Uttlt LV the P at ’ an, f ters used t0 ™ spectra. This expression * 
nas oeen fitted to a spectrum observed by the DMSP-F2 satellite (ref. 13), 
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THE CHARGE-STABI1.IZED POISSON ITERATION 
The Poisson-equation can be written as 
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potentials ind ?*’ ! ntl the charge density is determined using the 

potentials of the previous iteration. This method can be shown to be unstable 
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Equation (4) may be linearized about-the previous potential iterate 


q(^ v ) - q(tf v_1 ) + q'(^ v-1 ) * (tf v . ^v- 1 ) 


= a 9 /3 ^» and the ^ dependence has Deen dropped for claritv With 
this expression we may write the implicit Poisson iteration scheme 
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The finite element approximation to (5) produces the matrix equation 
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which is calculated to be-the maximum allowable charge for the element. 
Because of the artificial amplification argument, S is often the more real- 
istic total for an element. Before deriving S, we define the barometric 
potential i&t, « tn(n^) which is- the potential for which Q * 0 and note that 
it is important that S » Q as t> * Ab if quasi neutral regions are to be 
modeled correctly. To determine S, consider a capacitor with potential dif- 
ference (^5 - 4>) , area h% and a separation of h. The charge q c on this 
capacitor is given by 

"c * ClV - fi§~ (*b - *> kT 


In the units of our previous Q, q c becomes 
Q m ■ <*h(d b - 6) 


which is the maximum allowable charge per element, with the parameter a, 
adjusted to insure that is maximized.- . Thus at each node, we choose for 
the charge 

IS I - min.(|Q M |, |Q |) 

with 

/ -ah for S = Q m 

S ' = 

i -h\ - ^ exptf for S = Q 


The effect of this algorithm is this: If a problem has been specified 
where a boundary potential would be screened in less than a zone or two (the 
limit of any code's resolution), sufficient sheath charge will be redistrib- 
uted so as to allow the potential to be screened over the minimum number of 
zones that are consistent with stability. When this occurs it is necessary to 
have a modified criteria for locating the sheath "edge". Our choice is to 
place the edge at ^ = tn(ox z /h z ) (when > kT/e) which is the 
potential at which Q' n = Q'. This is the potential contour that marks the 
region where the most drastic charge reduction occurs. 

SHEATH PARTICLE TRACKING 

Internal to the sheath boundary, strong electric fields will cause sig- 
nificant bending of ion trajectories and focusing of currents. POLAR models 
these effects by tracking "particles" inward from the sheath edge to the 
object surface. This tracking provides both the distribution of surface cur- 
rents, i.e. currents to surfaces, and the ion density within the sheath. 

Prior to performing this calculation, we must know the suitable sheath 
location and the currents from the plasma to the sheath. 
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te^c&M ^r?^tta°l/?2! y ; the **® hdr ' Potential wrnJa 

fiai«r^ t inIl 4 l!^ t ?h t ! ,e *S ueStioa of sKeath boundary placement. The results of 
figure 8 indicate that the current density enhancement factors are not tremnn 

dously sensitive to the sheath boundary potentials for the rStwSrd SLles 
£jT H«t G ~ ward a °9 le !* bbere is a much greater enhancement sensitivity but this 
?K*iv Carr ^ ove r ,nt0 a sur> tace charging sensitivity. A lower Dotential 
sheath boundary will tend to extend ellipsoidally in the wake direction oro- 
ucing increased sheath surface area. Focusing within the sheath/wake area 
will collect these lower sheath currents and reduce the sensitivity to sheath 

J» IZlXITr positive c u rrln ? ’ sour c es 
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Before these particles are assigned currents, a test particle is ejXd out 

bw?iefrthlt wSold e shi»?H S * f f ?' ements “ P robe for objects or potential 
rill l hat 1? shield P^esheath currents from that plate. If such obsta- 
cles are formed, the particles from that plate are deleted. 

» Once all particles, and currents have been assigned., particles are 

fo?th C S?oJi°the t 7 e iiriJ raJecto r^® by a " pusherV This pusher sweeps back and 
forth along the Z axis, operating on successive single slices of elements ftwn 

sl c f/ potentials). Within each slice, trajectories advlice^n the X 
and Y directions until they reach the present Z or Z+l slice boundaries X Y 

Hn tbe V . ehiCle Sur f ac6 * If the P usher * sweeping i^the^z’direc- 
f x J tln ? at Z+1 are continued in the next push, whereas parti- 
.m C L fh 4 V 1 1 (moving in the -Z direction) are stored on disk and picked 
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RESULTS 

Two model calculations are presented, the quasisphere and the shuttle 
Fortoth problems, the plasma is an 0* plasma wiih kT . oa e!. i>. ^ 

speed * 0 ? ft T-fi'fem/c^i 80 ^ m ° de ] s . hav ® the objects moving at an ion Mach 
speed of 8(8 km/ sec)* Neutral ion density contour plots for these obiert* 

b® e o P r dsented in figures 4 and 5. All of the contour plots are 2-0 cuts 
nln?U 9h thil e 3 :° arra £ s of potentials or densities. In the quasisphere (Q-S) 
th f. CU J S 90 through the center of the object and its wake/ For the * 
shuttle, the longitudinal cuts (showing a side profile of the orbiterl run 
just to the side of center and through one of the engines. ' 

covered"wUh^ J he Q i U „" 10de,ed « ° grounded conductor 

covered with 0.1 mm thick kapton and a “quasiradius" of 2.5 m. At t * 0 the 

?ial * Thfo S U ?s°IlK% 7 ^ SWU J hed "• “I th a " sur fgces afzero ^,t - 
tial. The Q-S is allowed to charge for a total of 7.8 seconds in 13 intervals 
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of 0.6 sec each. Ambient electron spacec.harge factors and auroral electron 
surface currents are calculated at each step. All secondary electrons are as- 
sumed to escape. Ion surface currents are determined from sheath particle 
tracking at the 0.6, 2.4, 4.8 and 7.2 sec intervals, while the 0-0.6 sec in- 
terval uses the precharge estimates of ion currents. Figure 9 shows a poten- 
tial contour plot with- linearly spaced contours at 2.4 sec, where the- maximum 
surface potential is -130 volts on the woke side of the Q-S. Figui 10 is a 
contour plot of the same potentials us-ing logarithmically spaced c*.. :ours, 
where the sheath contour is chosen to be -0.6 volts and is indicated by the 
X's. A subset of the sheath ion trajectories are also shown. Figure 11 shows 
linear potential contours for the Q-S at 4.8 seconds with the new sheath loca- 
tion (-0.6 volts) indicated again by X's; the maximum surface potential is 
-252 volts. The end point of the calculation is shown in figure 12 at 7.8 
seconds where a rough equilibrium has been established and the maximum surface 
voltage is -447 volts. The apparent discharging of the side surfaces at 7.8 
seconds when compared with the Q-S at 4.8 seconds is a plotting illusion due 
to the same number of contours being spread over a large potential differ- 
ence. The long charging times reported here are not realistic, and are due to 
the combination of grounded conductor and thin dielectric which produces a 
high capacitance to ground. This fixed ground has been retained as a develop- 
mental convenience and will be allowed to float in the future. 

A preliminary model of shuttle orbiter charging is presented in figures 
4, 13, 14, and 15. The materials specified for this model are not realistic, 
and again, the ground was not allowed to float. Starting with the ion den- 
sities shown- in figure 4, the orbiter is exposed to the same plasma and aurora 
as was the Q-S with a mesh spacing of L.16 m/grid unit. The shuttle is al- 
lowed to charge for only 3 seconds in three steps. Potential contours at the 
3 second point are presented in figure 13 for the same cut as figure 5. The 
maximum surface potential is -101 volts. The lowest contour is at -20 volts 
so it is not possible to see lower space potentials in the wake, but one can 
observe the compression of contours on the underside compared to those in the 
cargo bay. Figure 14 is a cut through the ion density data at Z = 17 (.see 
figure 13). In figures 14 and 15, the complete projected silhouette is out- 
lined whereas only a portion of the cargo bay floor, wall and doors actually 
lie in this Cut. Since the plasma is flowing predominantly along the length 
of the orbiter, the ion density wake shown in figure 14 is due to portions of 
the orbiter upstream of the cut. Figure 15 shows potential contours in the 
cut at Z = 17, where significant charging can be observed on the outside Of 
the bay doors which lie in the wake of the wings, but not On the inside which 
was a grounded conductor. 

The shuttle model employed 41,175 grid points, but due to POLAR's disk 
data management and segmented construction, only about 70,000 words of core 
memory were required. On our UNIVAC 1100/80, the shuttle model calculation 
required about 3 hours of CPU time. 

CONCLUSIONS 

The results of these POLAR calculations are quite preliminary and are 
presented primarily to demonstrate the capabilities of the code. However, we 
are quite pleased with the close resemblance of the quasisphere model to the 
DMSP charging events reported in these same proceedings by Burke and Hardy. 
Perhaps the most significant result of these calculations is that POLAR's 


design criteria are proving to be realizable, and although more development 
will be required we anticipate that POLAR will become- a usef-wl scientific and 
engineering tool. 
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TABLE ~ POLAR /IONOSPHERIC. PARAMETERS 
Orbit Altitude 150-350 km 


Ambient Plasma 

Temperature (eVl 
Ion 0ensity(cnr 3 ) 

Ion Species 

Neutral Density (cm” 3 ) 


0. 1-0,3 

IQMO 5 + + 
0 * N .,He + , h + 
10 &- 10 11 


Ambient Electron Acoustic Speed 
Ambient G* Acoustic Speed 
Satellite Velocity 


V ae ■ /kl/m m 133 km/ sec 
V ao « 1 km/ sec 
8 km/ sec, M - V/V ao « 8 


Auroral Electrons 
Energies 
Fluxes 

Integrated Currents 


-10|-10| eV 
<10 s /cm‘*sec *str*KeV 
<100 iiA/nr 


Magnetic Parameters 


Magnetic Field Strength 

5 x 10” 5 Tesla 

(0.5 Gauss) 

Particle 

Gyroradius 

Gy r opened 

Ambient Electrons 

(0.1 eV) 2 cm\ 


Secondary Electrons 

(3 eV) 18cm ( 

7 x IO“?-sec 

Auroral Electrons 

(10 KeV) 4 m 


Ambient Ions (0 + , 0.1 eV) 
(H\ 0.1 eV 

/ 

3 m 
0.2 m 

2 x 10” 3 sec 
1.3 x lO” 3 set 

Cycloid Spacing Observed At 8 km/sec 
Electrons 0.6 cm 

H* 10 m 

0 + 160 m 



Characteristic Lengths 


Ambient Debye Length 

Xp < 1 cm 


Shuttle Orbiter 

L . 37 m 



Sheath Thickness For Potentials 1 kV 5 kV 

Around 3 m Radius Sphere* ** 6.9 m 11.6 m 

Planar Child-Langmuir + 9.32 m 31.2 m 

t Calculated for xp « 0.74 cm, kT « 0.1 eV 
** From Figure 72, Al'pert et ali (ref. 21) 
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•J-C a.0 7.0 12.0 17.0 22.0 27.0 32.0 37.0 42.0 <7.0 92.0-9.0 02.0 

z~ (grid units) 

Figure 3. Two-dimensional cross-section of POLAR's computational and object 
grids. The combination of orthogonal transformations of the object 
and object grid with staggering in the X*Y node slices allows for 
any orientation of_object and plasma flow. 
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Figure 4. Two-dimensional cross-section of normalized neutral ion densities, 
for a slice one unit from center, running through an engine. One 
mesh unit * 1.16 m, ion Mach speed * 8, contour interval * 0.11, 

M 1n(I _? ■ . . . ^ 





Figure 5. A central slice of ion denstt 
mesh uni ti= 1.0 m, ion Mach 
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Figure 


Figure 8 


differential electron flux 



. Fit and comparison of POLAR's analytic spectral form to an 
energetic electron event observed on DMSP-F2 (ref. 13). Assuming 
isotropy the total current carried by this spectra was 14 u A/nr. 



. Current density to a sphere of potential PHI * eV/kT in a Mach 8 
plasma using orbit limited theory, cos(e) « 1.0 for ram direction, 
cos(e) * -1.0 for wake direction. 
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Figure 9. 
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Figure 10. 
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Potentials and cross-section as figure 9 but with 

]ahf^Id hmiC 4 . C0nt0urs * The cor, tour interval is aloq(-V) - 0 87 
labeled contours are in volts. The lower eft tini/ibl 

the sheath location (-0.6 v), where a subsetn/l?^ I he . x s " 
illustrate the sheath ion tracking? ° f tra J ectori( * 
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Figure 11. 


z (meters.) 

Quasisphere potential contours after 4.8 secs of charging. Th 
maximum surface potential is -252 volts % and the contour inter- 
val is 28 volts. X's again mark the -0.6 volt sheath contour. 
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Figure 12. Quasisphere potential contours after 7.8 secs of charging, where 
some surfaces have reached equilibrium. The maximum surface poten 
tial is -447 volts and the contour interval is -49.7 volts. 
















